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FOREWORD 

This report was originally intended to be a United States contribution 
to the work of Study Group 6 (propagation in ionized media) of the 
International Radio Consultative Committee (CCIR) . The area of Extrater- 
restrial propagation is a new one to Study Group 6 (see Appendix A) so it 
was felt that a fairly substantial contribution was indicated. However, 

' ^n completed, the report was too long to be included in its entirety, 
id Section II was deemed too mathematical. Upon suggestion from Dr. Rush 
(chairman of U.S. Study Group 6) it was decided to reproduce the full report 
as a JPL Publication. 
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ABSTRACT 


The International Radio Consultative Committee (CCIR), in response 
to a U.S. initiative, approved a new Study Programme in 1978 entitled, 
"Radio propagation through the solar and other extra-terrestrial ionized 
nedia." The work described in the present report was carried out in 
response to the CCIR call for contributions in this area. It is intended 
for engineers with some background in radiowave propagation and frequency 
allocations. 

The report first treats the present S- and X-band communications 
needs in deep space to illustrate the aspects which are affected by 
propagation through extraterrestrial plasmas, and develops (in Table III) 
the magnitude, critical threshold, and frequency dependence of some 
eight propagation effects for an S-band propagation path passing within 
4 R of the Sun. 

The second section deals with the theory and observation of propaga- 
tion in extraterrestrial plasmas. The various plasma states along a near- 
solar propagation path are illustrated. Classical magneto- ionic theory 
(cold anisotropic plasma) is examined for its applicability to the path 
in question. Table IV summarizes the characteristics of the plasma states 
found along the path and indicates the errors in some of the standard 
approximations. 

The third section, on models of extraterrestrial plasmas, is primarily 
directed towards modelling the electron density in the solar corona and 
solar wind, but does also offer some cursory information on the terrestrial 
planets plus Jupiter. 
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SECTION I 


PROPAGATION EFFECTS IN DEEP SPACE COMMUNICATIONS 
THROUGH EXTRATERRESTRIAL PLASMAS 


1. INTRODUCTION 

The effects produced by Ionized media beyond the Earth's ionosphere on 
radio propagation are important in the selection of operating frequencies 
for future deep-space missions. The purposes of this report are twofold: 

(1) to identify those propagation parameters that are apt to be controlling, 
or at least significant, in the selection of frequency; (2) to examine those 
parameters in sufficient detail to assess the magnitude of their change and 
their behavior as a function of frequency. 

To accomplish this, the report contains a consideration of those 
telecommunications engineering and science requirements of space systems 
chat are most apt to be affected by propagation effects. The magnitude of 
the propagation effects for a sample case and their behavior as a function 
of frequency are then tabulated. Theoretical and observational aspects of 
propagation in plasmas are summarized in Section II and illustrated with 
examples of propagation through the solar corona and Interplanetary medium. 
Section III considers physical models of the solar corona and solar wind, 
and presents sample profiles of planetary ionospheres. 

One of the major concerns In the deep-space research program has to do 
with the use of yet higher frequencies (to the neighborhood of 30 GHz) for 
future deep-space missions. It raises questions both from the standpoint 
of the performance of the deep-space links considered in isolation, and 
from the standpoint of frequency allocations and sharing where the needs 
and characteristics of other services are involved. The Earth's atmosphere 
including hydrometeors will be a major factor, but its effect on deep-space is 
not significantly different from the geostationary satellite case. The effect 
of extraterrestrial ionized media, particularly in the solar corona, is unique 
to deep-space telecommunications. 

2. PROPAGATION EFFECTS RELEVANT TO DEEP-SPACE TELECOMMUNICATIONS 

The radio link connecting the Earth and the spacecraft has three 
primary functions. Telemetry is the first of these and has to do with the 
transmission of downlink data telling us about the spacecraft's surroundings: 
for example, video data of a planet's surface. Telemetry is characterized 
by its volume and the requirement for moderate to high quality data (e.g., a 
video picture contains enough redundancy that a few errors in data bit detec- 
tion will not be catastrophic, whereas loss of data on plasma structure would 
be much more serious). The second function is command, which refers to 
information sent to the spacecraft from the ground. This information con- 
tains directions to the spacecraft to alter its preprogrammed instructions 
and to take specified action with specified parameters at definite times. 
Command is characterized by low data volume and a requirement for extreme 
accuracy. The third function of the radio link is "radiometrics", a term 
that includes the development of information for navigation and scientific 
purposes from characteristics of the received radio signals. These charac- 
teristics are imposed by the position and velocity of the spacecraft, and 
by the propagation medium. 
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Two principal factors serve to alter the radio signals in their passage 
between Earth and spacecraft: the propagation medium, which may signifi- 

cantly degrade the signal for all three functions, and the velocity of the 
spacecraft relative to the ground stations. By exceedingly accurate 
measurements of the radiowave characteristics, the position of the spacecraft, 
its velocity, and its acceleration can be determined. Simultaneously, in- 
formation is acquired on the structure, composition, and temperature of the 
propagation medium (typically the solar corona or the atmosphere of a planet 
that may occult the signal, or the Interplanetary plasma state along the 
ray path), and on planetary and solar gravitational fields. These radio- 
metric measurements are characterized by the need for a high order of long- 
term frequency stability (for example, to parts in lO 1 ^ for the Earth station 
for gravitational waves), and by the extent of data processing required to 
turn the signal into information [Edelson et al, 1979]. 

The performance of these three radio-link functions is dependent on the 
signal-to-noise ratio (SNR), which may be defined by the equation: 


P r /N o " P t + G t + G r “ L 1 “ L 2 “ 20 108 * “ 10 log d - 10 log T g - 27.34 dB 


where 


P : transmitter power (dbW) 


( 1 ) 


G. : transmitting antenna gain with respect to an isotropic 

radiator (dBl) 


G r : receiving antenna gain (dBl) 

: attenuation in propagation medium different from free 

space (dB) 

L 2 : all other losses from receiving antenna to first stage of 

receiver (dB) 


f : frequency (GHz) 

d : distance in astronomical units (1 AU ■ 1.496 x 10 11 m) 


T 

s 


system noise temperature including additive noise (K) . 


An example illustrating the magnitude of these parameters for the 
case of the NASA Voyager spacecraft at the distance of Jupiter is given in 
Table I. The locations and diameters of the antennas of the three NASA Deep 
Space Network (DSN) complexes are shown in Table II. A comparison of the 
losses incurred simply due to distance, between satellites in geostationary 
orbit and spacecraft on planetary missions is shown in Figure 1. The 
nearest of the inner planets is seen to experience 60 dB more loss for any 
given frequency; the nearest ol the outer planets experience 87 dB more loss 
than for the geostationary case. 

Increasing the gain of a parabolic antenna (commonly used for deep-space 
work) decreases the beamwidth as illustrated in Figure 2 (where the effect 
of surface Irregularities is also portrayed). If there is angular broadening 
due to plasma effects this may spread the angular spectrum of Incoming plane 
waves to angles outside the 0.3 dB cone of the antenna thereby causing a 
decrease in signal level of 0.1 dB or greater. 
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The telecommunications parameters of a typical deep space mission are 
given in Table I and the antenna locations and dimensions in Table II. While 
the spacecraft transmitter power is seen to be modest, only 21.3 watts at 
X-band (8.4 GHz) and 6.6 watts (18.9 watts in the high power mode) at S-band 
(2.3 GHz), the input power required to generate it represents a substantial 
fraction of the total spacecraft power. 

Phase modulation is used to impress information on the carrier signal, 
and several components of the single-bit stream are separately coded to 
achieve the bit-error-rate required for each function or experiment. This 
scheme is called concatenated coding. In the case of Voyager, it allows 
a two-orders-of-magnitude BER difference between imaging and nonimaging 
data. The cost is in increased ground-processing complexity. Only recently 
have ground computer speeds increased to the point where decoding of 
concatenated coding becomes feasible. 

Criteria for determining just what constitutes a significant change in 
a propagation parameter are not easily established. A case in point might 
be broadening of the carrier signal that is used to synchronize the receiving 
system on the ground. The phase lock loop has a bandwidth of tens of Hertz 
(function of received power), so that broadening of the carrier by plasma 
fluctuations to say 50 Hertz could cause loss of lock, if the received power 
level were to be decreased simultaneously by amplitude scintillation. 

Similarly with the telemetry signal, if the margin has been exhausted, 
a further decrease of SNR by 0.7 dB will cause a ten-fold increase in imaging 
error rate, and 1000-fold for nonimaging error rates (assuming curves 3 and 
4 of Figure 3 for inner and outer coding, respectively). Error considerations 
in deep-space telecommunication systems are governed by the error source with 
the largest standard deviation. At the present time an uncertainty of 0.1 dB 
is significant. 

2.1 Radio Science 

In general, the effects produced by plasma on communications are 
undesirable. However, in radio-science investigations, these plasma effects 
while at times undesirable, as in the case of very precise celestial mechanics 
measurements are, at other times, used to develop information about the propaga- 
tion media. Both situations deserve comment. 
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TABLE I 


Typical NASA Voyager Telecommunications Parameters 


Parameter 

Unit 

Value 

Telemetry 

X-band S-band 

Command 

Carrier frequency, f 

Gi^iertz 

8.4 

2.3 

2.1 

Transmitter power, P 

dBW 

13.3 

8. 2 t 

40.0 

Antenna gain, G t 

dBi 

48.1 

35.1 

60.79 

Range at Jupiter, d 

AU 

4.55 

4.55 

4.55 

Receiving antenna 





diameter, D 

m 

64.0 

64.0 

3.7 

Receiving antenna 





effective gain, 

dBi 

70.2 

61.6 

33.7 

System losses, L^+L2 

dB 

0.8 

0.4 

0.0 

System noise 





temperature, T g 

K 

28.5 

22.3 

1500.0 

Power flux density at 

2 




receiving antenna, F 

dB(W/m ) 

-186.2 

-204.4 

-147.0 

Signal-to-noise ratio, 





P /N * 

dB(Hz) 

58.1 

43.6 

55.7 

r o 









, -3 

*115.2 kilobits per second 

requires a P 

./ N of 53.2 

dB for a BER 

of 5 x 10 

The 16 b/s command channel 

requires a Pj 

,/N° of 35.5 

dB for the required 

command bit error rate of 

5 x lO" 6 . 




tThe S-band transmitter power for the high-power mod 

e is 12.8 dBW. 
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TABLE II 


Antennas of the Deep Space Network (DSN) 


Location 

DSN No. 

Antenna Size 
(m) 

Frequencies* 

Canberra, Australia 

Tidblnbilla 

42 

I' 6 

S-band 

Tidbinbilla 

43 

64 

S-, X-band 

Honeysuckle Creek 

44 

26 

S-band 

Madrid, Spain 

Robledo 

61 

26 

S-band 

Robledo 

62 

64 

S-, X-band 

Cebreros 

63 

26 

S-band 

Goldstone, California 

Pioneer 

11 

26 

S-band 

Echo 

12 

34 

S-, X-bard 

Mars 

13 

64 

S- , X-band 


*S-band frequencies in deep-space use are 2100-2120 MHz uplink, and 
2270-2300 MHz downlink; 

X-band downlink frequencies are 8400-8440 MHz. 
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FIGURE 1 

Magnitude of basic transmission loss in free space for deep-space 
telecommunication distances. A * Geostationary distance, 

B ■ range of distances to inner planets, C ■ range of distances 
to outer planets. 
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FIGURE 3 

Bit error rate associated with different methods of encoding or adding 
redundancy to transmitted data plotted as a function of signal energy per 
bit (E^)» divided by noise spectral density (N^). (After Edelson ot al , 1979) 

(1) uncoded case, 

(2) code used in NASA missions 1969-1976. 

(3) NASA Voyager imaging code (inner code). 

(4) Voyager non-imaging science concatenated code (combined inner and 
outer codes). 

(5) Theoretical upper bound on performance Improvement (infinite hand- 
width case), Shannon’s limit. 



1,1 Celestial Mechanics 


The effects of gravity on a spacecraft orbLt causes detectable changes, 
which are used for harmonic analysis of the gravitation field of the planet. 
These analyses can then be interpreted in geophysical terms and related to 
the internal structure of the planet. 

On a more esoteric level gravitational waves are predicted by Einstein's 
general theory of relativity. If non-grav i tat ional wave fluctuations on 
the received frequency can be kept below 1 part In 10“' J , then It should be 
possible to detect gravitational waves of a 10-^-seeond duration or longer, 
assuming time standards of at least this accuracy are available. Either 
plasma effects on t he phase velocity need to he reduced below these levels, 
or a way needs to be found to determine the extent of the plasma contribution. 

General relativity also predicts a decrease in the group velocity of a 
radio wave due to gravity as it passes close te a massive object such as the 
Sun or .Jupiter. This Increase in apparent range is one incentive for improv- 
ing measurement of apparent range to decimeter or centimeter accuracies at 
distances of several All. 

1.2 Plasma Measurements 

Two types of plasma information are derived from the telooommuni- 
cat ion (or special beacon) signals. The first is from the spectral 
broadening, polarization change, phase and amplitude scintillations, 
and angular broadening. These measurements are all concerned with 
fluctuations in the plasma. The second, having to do with bulk prop- 
erties includes informal ion on total electron content which is derived 
from range and doppler data and Ear ad ay rotation measurements which 
have been used to gain information on magnetic fields. The principal 
targets for these measurements are the solar corona, the solar wind, 
and planetary ionospheres . 

Spectral broadening and phase scintillations are two of the most useful 
radii’ scattering phenomena that can he observed with monoc hromat ic , coherent , 
and point source spacecraft signals. Spectral broadening is only useful near 
the Sun wheie the induced broadening exceeds the inherent line width of the 
radio system, which is typically 0.2 Hz. Since spectral broadening measures 
t requeue les higher than 0,2 H/., it provides information only on the sma Mot- 
scales of the inhomogene i t v spectrum las also do amplitude scintillations, 
e.g., .lokipti [l‘)7l]l. l’hase idopplorl sc i nt i l 1 at ions , however, provide 
informat ion on the full range of scale si.es. 

Section II contains a discussion of these eftects; also, a good and 
readable descript ion of those topics is found in Woo |l l )7/|. 


2 Telemetrv 


The most critical telemetrv parameter in practice Is s i gna 1 -t o-no i se 
ratio, and tor that reason it is worth examining, in some detail. Shown in 
Figure ' are the 1’ i t error rates associated with ditferent methods of 
encoding data where the abscissa scale is energy per bit U'g) divided by 
noise spectra! densitv (N^,) . The energy per bit is t lie received power l’ 
t i nuns the h i t duration T^. In the case of a modern mission such as NASA’s 

Vovagor, the error sensitive data is protected bv an outer code betoro 
being combined with error- i nseus i l i ve data and put through the intur 
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encoder. The resultant is satisfactorily robust, and mission requirements 
are satisfied by E /N * 2 or 3 dB. However, the bit error rate decays 
rapidly on both curves 3 and 4 of Figure 3 and this situation changes from 
sublime to intolerable if the SNR decreases by 2 dB. 

Good engineering practice dictates the provision of a power margin. In 
conventional satellite systems the power margin represents a reserve 
capacity that is not available to space systems which are stretched to their 
observing limit. Therefore, there is an argument for not considering the 
margin in determining the effect of a decrease in signal level due to prop- 
agation factors. 

Let the criterion for serious deterioration of signal quality be a 
decade increase in bit error rate at the limiting threshold. For imaging 
telemetry this threshold is at a BER of 10 -3 . From Figure 3 curve 3 a 
change from a BER of 3 x 10 “* to 3 x 10“ 3 is caused by 0.66 dB decrease in 
E r /N . Similarly, for nonimaging telemetry (curve 4), the corresponding 
decade change in BER around 10“ 3 is caused by a 0.19 dB decrease in Eg/N Q . 

2.3 Navigation 

Interplanetary navigation is consistent with traditional navigation in 
that it involves a map, a travel plan, a means of periodically determining 
the ship's location and course, and a method of selecting a new route when 
the ship has drifted off course. The map of the solar system is constantly 
being refined. Missions to near and remote planets are planned and space- 
craft trajectories determined. The location and velocity (six numbers) of 
the spacecraft are determined en route (by radio means), and, as a result, 
corrections are made that allow the spacecraft to reach its interplanetary 
goal. 

The map of the solar system has been materially improved by planetary 
radar measurements. As late as 1960, the astronomical unit, the mean Sun- 
Earth distance, was in error by 66,000 km. The present value of 149,597,871 
km is known with an uncertainty of about 1 km [Melbourne, 1 976 j . Radar 
range measurements of the terrestrial planets are accurate to about 50 meter 
under favorable conditions. 

Two methodo fnr establishing the position and velocity of a spacecraft 
have been used long enough to be termed classical: the measurement of the 

round-trip travel time of a radio signal (range), and the measurement of 
the Doppler shift of the radio wave after it makes the round-trip to the 
spacecraft and back to Earth (rate of change of range, or radial velocity). 

2.3.1 Range 

The ranging system presently employed for NASA deep-space missions 
switches the phase of an S-band carrier back and forth by a radian, say, 
every microsecond. This is equivalent to a 0.5 MHz square-wave phase 
modulation. A microsecond corresponds to 300 m, one way, or 150 m range; 
and it is necessary to decide which 150 meters contains the spacecraft 
(i.e., there is an ambiguity modulo of 150 m) . This is done by successively 
halving the modulating frequency up to 20 times until the ambiguity is 
clearly resolved. 
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Range accuracy by the "classical" method has been about 3 meters In 
one way range [Melbourne, 1976}. Through careful consideration of the 
various error factors 1 to 2 meter range accuracy has been achieved at 
S-band in the present NASA Deep Space Network (DSN) . Further improvement 
in S-band ranging for the inner planets is limited by solar plasma scintilla- 
tion [Anderson and Estabrook, 1979]. 

A long-range plan is to move to X-band for ranging. The expectation 
is that 50 cm ranging at Jupiter distance will be available in ten years, 
and 10 cm ranging at a later date. This improvement is based on the wider 
bandwidth available at X-band and the expected f~^ decrease in plasma effects 
at 8.4 GHz as compared to 2.2 GHz. 

2.3.2 Range change 

A radio signal of frequency f c sent to a spacecraft from the ground, 
turned around by the spacecraft transponder, and retransmitted to Earth 
would arrive with a frequency f^ which is given by f Q - f; ■ (V r /c)2fo, 
(typically f 0 - f-| ^ 100 kHz), where V r is the radial velocity of the 
spacecraft measured in the outward direction; in this example, it is 
7.5 km/s. If the frequency can be determined to 0.01 Hz, velocity can be 
determined to 0.75 mm/s. The actual measurement is in terms of phase and 
the accuracy depends on the integration time. 

These two values of range and range change are quite adequate for 
navigation in one dimension. A problem arises with respect to the other 
two coordinates in the plane orthogonal to the radial direction, namely 
declination and right ascension on the celestial sphere. Here the accuracy 
of range is not much help. For example, with a spacecraft at a distance 
of 0.4 AU, a displacement of 600 km perpendicular to the Earth-spacecraft 
direction results in an increase of only 3 meters in range. As a conse- 
quence, a different family of techniques has been developed to determine 
position and velocity in the plane perpendicular to the radial direction. 

2.3.3 Angle determination 

Differenced range is a method employed in the NASA Deep-Space Network. 
Its features are illustrated in Figure 4 where A p = B sin 6 is the differen- 
tial range. As shown in Table II DSN complexes are located in California, 
Australia and Spain. The U.S. and Australian complexes are located at 
approximately equal north and south latitudes but differ in longitude by 
about 90°. Hence the angle measured by differenced range has a large right 
ascension component as well as a large declination component. As only one 
tracking station can be transmitting to a spacecraft at a time, simultaneity 
is actually impossible. Instead the range measurements are "nearly simul- 
taneous." For example in the case of a Saturn conjunction Goldstone, 
California and Canberra, Australia have the spacecraft in both antenna beams 
for 4-1/2 hours. The round trip transit time for a radio wave to Saturn 
and back is 3 hours leaving 1-1/2 hours when either station can command the 
spacecraft and still receive a response back. During that period command is 
handed back and forth several times from Goldstone to Canberra in order to 
achieve near simultaneity. 
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FIGURE A 


Differenced range as an argle measurement. 


2. 3. 3.1 Very long baseline Interferometry (VLBI) 

>. new technique for angle measurement Is being tested In the Voyager 
mission. This Involves very long baseline Interferometry (VLBI), a 
technique developed by radio astronomers [Cohen 1969]. A spacecraft signal 
modulated by a pure sinusoidal tone, or an on-board noise source Is 
observed by two DSN stations with synchronized clocks. Accuracy of frac- 
tions of seconds of arc can be achieved In short times. The allied technique 
of AVLBI finds the spacecraft angular position with respect to cosmic radio 
sources. Accuracies of about 10~^ arc seconds should be available from 
VLBI In the 1980s, and AVLBI promises even greater accuracies [Anderson and 
Estabrook, 1979]. 

2.4 Conmand 

The command link (Earth-to-space) is characterized by the need for 
extreme reliability (BER _< 10"5). Commands must get through despite pointing 
errors on the spacecraft, and weather effects on Earth. As a result, the 
link budget is established for almost omnidirectional receiving antennas on 
the spacecraft, and the command frequency is firmly established at S-band 
(in order to stay below the frequencies that become vulnerable to terrestrial 
weather effects). The data rate requirements for the command link are 
determined by the size of the onboard computers, which control the spacecraft. 

Host onboard activities are preprogrammed into these computers by 
periodic command sequences. It is desirable to complete these transmissions 
within the view period of a single ground transmitter (about 8 hours). To 
have the capability of completely reloading the computer in this time, the 
command link operates at 16 bits per second in the case of Voyager 
[Edelson et al, 1979]. The command signal (S-band) will be the last to be 
impacted by plasma effects, but it will be impacted when, fcr example, the 
spacecraft passes behind the solar corona at which point it suffers more 
degradation than the X-band signal due to the greater plasma effects at the 
lower frequency. 

A typical command signal consists of Manchester-coded data bits PSK 
modulated onto a square wave subcarrier. Manchester coding is used to pro- 
vide bit synchronization. The command channel error rate follows curve 1 
of Figure 3. 

2.5 Noise Considerations 

Natural noise at microwave frequencies is treated from the standpoint 
of emission noise in the Earth's atmosphere in CCIR Rep. 683 and Rep. 720. 
Noise temperature in deep space receiving equipment (including antenna 
losses) is treated in Rep. 536-1. System noise temperature for Earth 
stations at S-band runs about 22. 3K (see Table I), which at X-band rises to 
about 28, 5K. These extraordinarily low temperatures make the deep-space 
system particularly vulnerable to gaseous or plasma emission noise. 
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3. SUMMARY 


The telecommunications signal in deep-space missions serves a number 
of distinct purposes. Each of these uses has its own criteria for what 
constitutes significant signal degradation. Examples for the Voyager mission 
follow. 

Uplink 

The command signal increases its bit-error-rate (BER) a factor 
of 10 with a 1.1 dB decrease in signal-to-noise ratio. 

Downlink 

Radio science celestial mechanics experiments need to determine 
range to decimeters, and hopefully, centimeters at planetary 
distances. 

The information transmitted by the telemetry signal is degraded 
by a factor of 10 in BER by a 0.2 dB decrease in SNR for the 
non- imaging data and by 0.7 dB decrease in SNR in the imaging 
data. 

The navigation determinations currently rely on a minimum of 
3 meter accuracy in the doppler and ranging determinations and 
anticipate further accuracies of less than 1 meter. 

The estimated effects of ionization in the solar corona and solar wind 
for a spacecraft on the other side of the Sun, whose propagation path comes 
within four solar radii (4 R 8 ) of the center of the Sun, are tabulated in 
Table 111. It is clear that several effects have become of significant 
magnitude. Also note that not all the effects are inversely dependent on 
the frequency of the source. Those due to scattering can be significantly 
shallower in their decrease with increasing carrier frequency. 
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TABLE III 




Magnitude, critical threshold and frequency dependence at S-band of the 
various propagation effects for a propagation path traversing the solar 
corona to a closest approach of 4R, of the Sun. One way traversal. 


Effect 

Magnitude 
(at 2.3 GHz, 4R # ) 

Critical 

Threshold 

Frequency Dependence 
(f) 

(see Section II) 

Faraday Rotation 

20-200° 

(3-7° in ionosphere) 

0.5“ 

(see note 1) 

1/f 2 

Group Delay 

17 ysec 
(250 nsec in 
ionosphere) 

20 ns 

(see note 2) 

1/f 2 

Angular Broadening 

0.02-2 min arc 

+2 min of arc 
(see note 3) 

1/f 2,2 to 1 /f 3 

Spectral Broadening 

7-16 Hz 

0.2 Hz 
(see note 4) 

1 /f 1 * 2 to 1/f 2 

Absorption 

negligible 

0.1 dB 

1/f 2 

Dispersion 

10 nsec /MHz 

2 ns/MHz 
(see note 5) 

1/f 3 

Intensity 

Scintillation 

saturates (m ■ 1 ) 

m - 0.0024 
(see note 6) 

1/f 1 * 42 to 1/f 1,2 

Phase Scintillation 

(see 3.1 of 
Section II) 

+4° 

(see note 7) 

1/f 1 


Notes: 


1. Telecommunication signals are normally circularly polarized. Faraday 
rotation is used as a sensing tool for studying the solar corona. 
Measurement accuracy is to be better than 0.5° for effective system 
time constants 30-1 00s [Volland et al, 1977]. 

2. An accuracy of 3 m in range requires an accuracy of 20 ns in group 
delay. This is achieved in practice through careful combination of 
several techniques and correction for plasma effects [Melbourne 1 976 j „ 

3. This value (2 min of arc) is for the 64 m antenna at S-band and 
represents 0.1 dB decrease in signal level if the energy is spread 
uniformly over a +2' cone of angles. At X-band the corresponding 
value is +5 min cf arc. 
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Notes (continued) 


4. The value of 0.2 Hz is a lower limit on reported spectral broadening 
(Helios 1975 as reported in Woo 1978). In principle spectral broad- 
ening by a few tens of Hertz could cause the phase lock loop to lose 
lock but in practice this is unlikely. 

5. At present a 10 MHz band is available for ranging. A value of 

2 nanoseconds per megahertz corresponds to a difference in group delay 
of 20 nanoseconds across this band which in turn corresponds to a 

3 meter ranging error. 

6. A variation in power of +2.4% equates to +0.1 dB. The scintillation 
index is defined as the RMS fluctuations of intensity about the mean, 
relative to the mean. The +2.4% has been equated to this RMS 
fluctuation. 

7. Phase is accumulated with an accuracy of 4° by che Doppler-cycle 
counter in the range-change navigation system [Melbourne, 1976]. A 
more serious effect has to do with the ability of the phase lock loop 
to track rapid phase excursions. 
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SECTION II 




THEORY AND OBSERVATION OF PROPAGATION IN EXTRATERRESTRIAL PLASMAS 


1 . INTRODUCTION 

Plasmas occur in various forms and it is ne'essary to specify several 
parameters to determine the regime of the plasma. These parameters are 
temperature, magnetic field, degree of ionization, electron density, and 
homogeneity. The effect of the plasma in each regime on an electromagnetic 
wave passing through it will depend on the frequency of the wave relative 
to the local plasma frequency and the gyrofrequencies of the charged particles, 
the angle of the direction of propagation relative to the magnetic field 
(if any), and the polarization of the wave. A description of radio waves 
in plasmas is found in Report 222-4, and fuller descriptions may be found 
in Ginzburg f 1 9671 , Stix [1962], Shkarofsky et al. [1966], Wait L 1 968] and 
Yeh and Liu [1972]. 

The distinction between a warm and cold plasma is whether the wave- 
length iu the plasma is much greater than (cold) or much less than (warm) 
the Debye length, X^. 


/T\^ ^ 

X D “ 69 \N/ meters 

where T is the temperature in kelvin and N is the electron density (m ). 

In the Earth's ionosphere, the Debye length is typically 0.3 cm in the day- 
time and 1 cm at night, and the plasma may be considered cold at S- and X- 
(but not K-) band. In the solar corona, the Debye length is typically 0.1 
to 1 meter, and the plasma appears lukewarm to the 13-cm waves used at 
S-band for deep-space telecommunications, and warm for X-band and above. 

The critical magnetic field, B c (that is, the level at which a plasma 
ceases to behave isotropically), may be taken as the field strength, H, at 
which the energy in the magnetic field (1/2 y 0 H^) is equal to the thermal 
energy of the electrons (3/2 N k T) . Hence, 


B c • J3 N k T " 7.21 x 10 -11 ✓nT gauss 

—23 

where k is Boltzmann's const. (1.38 x 10 J/K) and y is the free-space 
permeability. ° 

In the Earth's ionosphere, B is typically 2 to 3 x 10 gauss, whereas 
the Earth's magnetic field at ionospheric heights is 0.3 to 0.6 gauss. 

Hence, the Earth's ionosphere is strongly magnetically controlled. In the 
solar corona, B c is ^lO - * gauss at 4 solar radii from the Sun's center (4R ) 
and the static magnetic field is -*10“^ gauss, so magnetic control is still® 
the rule. Illustrated in Figure 5 are the plasma domains traversed by a 
radio signal from a spacecraft near solar conjunction received at the 
Earth's surface (see Figure 6). We see that the 2.3 GHz signal is largely 
in ttie "warm-anisotropic plasma" domain until it passes into the "cold- 
anisotropic" domain (i.e., the Earth's ionosphere). This change of domains 
indicates caution in use of the cold plasma theory. Fortunately, the warm 
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FIGURE 5 

Plasma domains traversed by 2.3 GHz spacecraft, signal. Points A through E 
are specified in Table IV. Spacecraft located as shown in Figure 6. Domains 
are determined by the Debye shielding distance and the magnetic critical 
value B c (see text) and are; I - cold isotropic, II - cold anisotropic, 

III - warm anisotropic, and IV - warm isotropic. Plotted points are those 
marked on Figure 6: A - spacecraft location, B - solar corona (4 Rq) , C - 

Earths' ionosphere (300 km), D - plasmasphere (4 R_) , E - geostationary 
orbit (6.6 R_). W 



SPACECRAFT 



FIGURE 6 

Propagation path used for sample calculations. Spacecraft is at A, 
closest approach to Sun is at B, peak of Earth's ionosphere is at E 
Also shown are points D, just within the plasmapause, and C, the 
geostationary distance. 




v*. cold plasma distinction becomes import f it when the wave velocity of the 
signal becomes comparable to the particle velocities in the medium [Yeh 
and Liu, 1972], and that is not the case for any of the microwave links 
under consideration. 

The effect of magnetic control will be dealt with as the need arises 
in the topics treated below. 

In. the discussion that follows the following homogeneous plasma effects 
on radio waves are treated first: group delay, absorption, dispersion, 

Faraday rotation and Doppler shift. This is followed by a discussion of 
effects due to lnhomogeneitles , including phase scintillations, intensity 
scintillations, spectral and angular broadening, and coherence bandwidth. 

The propagation path used as an exampje in the following discussions 
is shown in Figure 6. It is the same one used in connection with Table III 
in the main body of this report. Table IV contains a listing of th^ 
various parameters used in the text for three locations, which are matked on 
Figure 6. An objective of Table IV and the discussion to follow it to 
illustrate circumstances where it is safe to use simple approximations. 

MICROWAVE PROPAGATION IN A HOMOGENEOUS ANISOTROPIC PLASMA 

The Appleton-Hartree equation for the complex refraction index in a 
homogeneous anisotropic ionized medium is giver, by 


n" - 1 - 


1 - 1Z 




1 2 
- Y /(I 
2 T /l 


tz,] ± jp 


/ 4(1 - X - IZ) 


FT 


( 2 ) 


where the nomenclature is that given by Ratcliffe [l959]j 



e,m charge and mass of the electron 

0 Angle between the direction of phase propagation and the external 
magnetic field 

v Collision frequency between electrons and other particles. 

2 2 
u> The square of the angular plasma frequency of the medium (Ne /i: Q m), 

where N is the electron density, and c ( „ is the permittivity of 

free space. 

f N Plasma frequency (w /2n). 

ui„ Angular gyrof requency or cyclotron frequency of the electrons 

(M i 0 / n) - 

w Angular frequency of the wave 

X (w N / w ) 2 . 

Y [e/ (raw) ]3 q . 

Y (Djj/w. 

Y Y cos 6 

y t Y sin e 
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Magnetic control, B/B | 0.66 (weak) 4.6 (moderate) 43.5 (strong) 40 (strong) 160 (very strong) 




R Wave polarization 


The implications of this equation (for a cold plasma) and the 
associated theory can be found in many excellent texts, for example, 
Ratcliffe [1959], Pudden [1961] and Davies [1968]. Its limitations and 
extensions may be found Wait [ 1 968 3 in Yeh and Liu [ 1 972] and Booker [1975]. 

Longitudinal propagation (0 = 0°) and transverse propagation (0 = 90°) 
offer significant simplifications, and have been worked out in detail. The 
quasi-longitudinal or Q.L. approximation [Booker, 1934] is possible when 

| yJ/2Y l | « | 1 - X - iZ | (3) 

This leads to the Q.L. approximation of (1): 


2 


n 


X 

1 - iZ + Y ’ 

jL 


R = + i 


(4) 


where two counterrotating circularly polarized waves are allowed. As the 
left side of (3) varies as fjj/f (and fjj, the electron gyrof requency , is 
normally 1 MHz or less), and the right side of (3) approaches unity with 
increasing frequency, it is clear that for frequencies above 1 GHz, the 
quasi-longitudinal approximation is good except for a small fraction of a 
degree around 0 ~ 90°. In the discussion that follows, the Q.L. approxi- 
mation will be assumed to be valid. 

2.1 Group Delay 

Group delay is discussed in Section 4.2 of CCIR Report 263-4 (and in 
the texts referred to above) for the isotropic collision-free case. It is 
instructive to examine the radio science and navigation needs in deep-space 
research to see if that discussion will suffice. 

The group refractive index iig may be defined in terms of the phase 
refractive index n 


n 


g 


d_ 

do> 


(nw) = n + w 


dn 

doi 


(5) 


This expression is evaluated for the collisionless case in Davies 
[ 1 968 J , and for the complete Appleton-Hartree expression (2) in Budden 
[1961]. 

If the Q.L. approximation (4) for the phase refractive index is 
reduced to the collisionless isotropic case (n^ ■ 1 -X), then from (5) 


n 

g 



n 


(isotropic collisionless case) 


( 6 ) 




which is a familiar case in which nn g - 1. However, retention of the O.L. 
approximation by substituting (A) in (6) and performing the differentiation 
in (6) gives 


n 


g 


1 

n 


i+ x(z t Y) - , 

2(1 - Z - Y) 


(Q.L. approximation with 
collisions) 


(7) 


Expression (6) is shown numerically in line 13 and expression (7) is 
evaluated in line 14 of Table IV. The difference between the isotropic 
value of the refractive index, n 0 , and that using Q.L. approximation 
is seen to be negligible, as is the error in computing n g , the group 
refractive index, as 1/n (the reciprocal of the phase refractive index as 
in (6)). Hence, for purposes of computing group delay, the approach used 
in Report 263-A is still valid for deep-space problems. 

The time delay, T (in seconds), of a radio signal traversing the deep- 
space-to-Earth path L (spacecraft to Earth station, shown in Figure 6) is 
given from (6) by 



_3 

where N is electron density in m , and f is the frequency in Hertz which 
leads to: 


4I . d , 

where AT is the group delay due to the presence of free electrons measured 
in seconds. It is a function only of the integrated electron density along 
the propagation path and the radio-wave frequency. 

As seen in Figure 7, the excess time delay for an S-band uplink and 
X-band downlink reached 12 ps. Assuming identical propagation paths, the 
delay on the downlink (8. A GHz) would be 7.33% that of the S-band (2.3 GHz) 
uplink. Thus, 11.2 ps can be attributed to the S-band uplink. From (9) 
we have 


I: 


TEC - I N dt * 7.493 x 10 6 AT f 2 el m 2 


= 4.4A x 10 hexem 


( 10 ) 


where TEC stands for total electron content, and the hexem is a unit of 
content defined as 1 hexeni * 10^ electrons/meter . 

The simple approximations of (6) through (10) will suffice for all 
conceivable microwave applications. 
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TIME SCALE 


FIGURE 7 

Range delay experienced by Viking spacecraft during January 1979 solar con- 
junction for S-band uplink and X-band downlink. Elongation is 1.02° over 
southern heliographic pole (R % 4 Rq per figure 6). Model based on 1976 
Viking solar conjunction. (Courtesy J. C. Breidenthal and T. A. Komarek, 


| 

| 




2.2 Dispersion 


A dispersive medium is one in which the refractive index (n) is a 
function of frequency (or free-space wavelength). In a plasma where 
to >> ujjj, and u> >> the dispersion is said to be normal if the phase 

refractive index is an increasing function of frequency. Differentiating 
(9) yields a particularly useful relation (assuming the path does not 
change) : 


dT 

df 


2.689 x 10 


-7 



seconds /Hz 


(ID 


For example, using the TEC value derived in (10), we find dT/df = 

9.813 x 10“ 15 s/Hz or 9.813 nanoseconds/ifflz for the difference in delay 
time for components of the signal spaced one megahertz apart in our example 
(approximately the configuration of Figure 6). 

2.3 Absorption 

It is known that a negligible role is played by plasma absorption in 
microwave propagation through the solar corona and interplanetary medium, 
but it is worth reviewing why absorption Is negligible. The path absorp- 
tion, L a , is given [Davies, 1968] as: 


L 

a 



decibels 


( 12 ) 


where A = 1.15 x 10 1 Nvf “ decibels/km (13) 

-3 

s is distance in km, N is electron density in m , v is the pertinent 
collision frequency, and f is the radio wave frequency in Hertz. The 
pertinent collision frequency turns out to be for coulomb collisions of 
electrons with ions [Rishbeth and Garriott, 1969]: 

V ei = ( 59 + 4 ‘ 18 1 o 8 10 (t3/N) } 10 ~ 6 nt ” 3,/2 seconds" 1 (14) 

-3/2 

Here the important point is that v e ^ depends on T . In other words, the 
collision frequency decreases with increasing temperature. As a conse- 
quence, at point B of Table IV (the closest approach to the Sun for the 
propagation path of Figure 6), the absorption coefficient from (13) at 

2.3 GHz is only: 

A = 3.2 x 10 11 decibels/km 
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If this absorption coefficient were to obtain for the entire propagation 
path of two astronomical units (300 million kilometers), the path absorption, 
L a , would amount to only 0.01 dB. Alternatively, if the absorption coeffi- 
cient were to increase by a factor of one thousand (to 2.5 x 10'-* m”^) as 
would be appropriate to 1.4 solar radii rom the center of the Sun, and if 
this absorption coefficient were to hold for 300,000 km, the path absorption 
would still only amount to 0.01 dB. Clearly, absorption at microwaves is 
negligible for the solar corona and solar wind. 

2.4 Faraday Rotation 

The two slightly different refractive indices given by (3) for the 
quasi-longitudinal approximation indicate that a linearly polarized wave may 
be treated as two counterrotating circularly-polarized waves with slightly 
differing phase velocities. For microwave frequencies, the binomial 
expansion may be used on (3): 

n = J 1 - X/(1 + Y l ) =1-1/2 X(1 + Y l ) 

which leads to the Faraday rotation formula for rate of change of the angle 
of rotation of the plane of polarization with distance: 


g.lxY L - 1.3548 x 1°- 10 NB l d eg /» 


(15) 


where N is electron density (m - ^), and Bl is the longitudinal component of 
the magnetic field in gauss. Good discussions of Faraday rotation measure- 
ments to determine total electron content in the Earth’s ionosphere are 
found in Garriott, et al. [1970] and Budden [l96l]. 

Faraday rotation measurements of radio paths chrough the solar corona 
are treated by Levy et al. [1969], Stelzried et al. [1970] and Volland 
et al. [1977]. 

The magnetic field emanating from the sun is roughly radial at distances 
greater than about 2 solar radii with a spiral effect dominating at longer 
distances. Parker [ 1 963] has approximated the angle of the spiraling lines 
by 


tan a 



(16) 


where r = distance from center of the sun, km 
fi s = Sun rotation rate, rad. sec" 
v - velocity of the solar wind, km sec 


As a consequence the magnetic field is about 1.7° off radial at 4Rg for 
the parameters given for location B in Table IV and about 46° off radial at 
1AU for the parameters for location C, so there is a slight asymmetry to the 
magnetic field as encountered by a radio wave following a path such as 
illustrated in Figure 6. If the solar magnetic field were perfectly radial 
(and maintained the same polarity) and the electron density spherically 
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symmetric, Faraday rotation would accumulate (wind up) in one sense of 
rotation up to point B of Figure 6 and then unwind by an identical angle 
and no rotation would be observed on Earth. This, however, is not the 
case as illustrated for Helios in Figure 8. Actually the magnetic field 
will have segments ol constant polarity with neighboring regions of opposite 
polarity thus greatly increasing the chances for non-cancellation of 
integrated Faraday rotation. Observed Faraday rotation for solar conjunc- 
tion of spacecraft has been analyzed and found to be in good agreement with 
current understanding of magnetic sectors in the corona and solar wind 
[Volland et al., 1977]. 

2.5 Doppler Shift 

The Doppler shift, Af, of a signal transmitted with a frequency, f, 
from a moving spacecraft (Figure 6) is given by 


Af 


f dP = 

. ii. 

c dt 

c dt 1 


ds 


(17) 


where P is the phase path from transmitter to receiver and n is the phase 
refractive index. As indicated in lines 11 and 12 of Table IV, the differ- 
ence between the isotropic and Q.L. approximation is negligible, and the 
phase refractive index n may be approximated by n 2 1 - X/2 just as can the 
group index n„ = 1 + X/2. As range change is acquired in tracking both 
from integrating the Doppler observations (D4> ) and differencing group path 
measurements (Dg), the difference (Dg - D<j>) is a measure of change in 
integrated electron content along the ray path [Callahan, J975]. This 
differenced range versus integrated doppler (DRVID) is used by NASA to 
provide one measure of the change in total electron content (TEC) along 
the propagation path. The phase path data utilized in DRVID are also the 
source of much of the data for the fluctuation studies discussed in the 
next section. 


3. MICROWAVE PROPAGATION IN INHOMOGENEOUS INTERPLANETARY PLASMA 

Interplanetary scintillations (IPS) have been studied using natural 
stellar sources since the discovery of IPS [Hewish, et al., 1964], The 
observing frequencies are commonly in the VHF and UHF band [Coles, et al., 
1974], but there are compensating advantages to using spacecraft signals. 
For example, in addition to the study of intensity scintillations, one may 
observe spectral broadening and phase scintillations if one has access to 
coherent monochromatic spacecraft signals. This aids considerably to the 
range of scale sizes of the medium that can be studied as shown in 

Figure 9 [Woo, 1977], The plot is in terms of the spatial wave number, K, 

which is related to the scale size l by K = 2n/£. Hence it can be seen 

that intensity scintillations are limited in applicability to the stuay 
of scale sizes from 60 m to 600 km, whereas phase scintillations extend 
from 1 km to several AU. 
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FIGURE 8 

Polarization angle (coronal Faraday rotation) measured on day number 241 of 
1975 during the second of two Helios-1 occultations. Recording stations 
are Goldstone, California and Effelsberg near Bonn. F.R.G. (after Volland 
et. al. 1977). 

: 
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FIGURE 9 

Range of heliocentric distances and spatial wavenumbers, K (related to the 
scale size H by K * 2 it/JI) for spacecraft radio scattering observations for 
2.3 GHz signal and SNR of 40 dB. P * phase scintillations, I ■ intensity 
scintillations, B ■ spectral broadening (after Woo, 1977). 




All aspects of scintillation become critical as the propagation path 
passes near the Sun and these will be considered in turn. 

3.1 Phase Scintillations 


Phase scintillations of spacecraft signals have been used to study the 
solar wind and solar corona because they contain information on a wide 
range of scale sizes [Woo, et al., 1976a, b; Woo, 1977; Woo and Armstrong, 
1979]. Extensive phase difference scintillation observations (2700 hours) 
have been made simultaneously at 2.3 and 8.4 GHz in 1976-1978 for helio- 
centric distance ranges, R, of 2-215 R # . (The symbol R # stands for the 
solar radius, a distance of 696,000 km, and 215 R^ is 1 AU, while R is the 
closest approach to the Sun of the propagation path and is sometimes called 
the impact parameter.) 

The 8.4-GHz signal is phase coherent and exactly 11/3 the frequency 
of the 2.3-GHz signal. Hence by recording the phases <J>i - 3/11 $2» where 
$1 and $2 are c ^ e phases for 2.3 GHz and 8.4 GHz respectively, the effect 
of phase change on the uplink signal is cancelled out. An example of this 
phase difference measurement is shown in Figure 10 exhibiting both small- 
scale (minutes) and large-scale (hours) variations. 

The phase difference spectra W^<j(f) are computed by Fourier transforming 
the time series, squaring, and smoothing the result [Woo, 1975; Woo, et al., 
1976a]. The resulting spectrum is shown in Figure 11. 

The orbiters of the Viking Mission to Mars were equipped with S- and 
X-band dual-frequency (2.3 and 8.4 GHz) coherent radio systems [Michael, 
et al., 1977]. Spectra for various elongations (Sun enter-Earth Station- 
spacecraft angles), c (one solar radius corresponds to an elongation of 
about 0,27 degrees), are shown in Figv 12 [Woo and Armstrong, 1979]., The 
high frequency instrumental phase noise is white and is about 10“^ rad^/Hz 
for the closed-loop spectra (inner three) of Figure 10 but is about 
10" 4 rad^/Hz for the open-loop spectra (outer three) where it had become 
necessary to open up the receiver bandwidth due to spectral broadening. 

A frequency spectrum of electron density, made up of in situ measure- 
ments made by a spacecraft as it moves along its orbit, is a one-dimens ionaJ 
wave number spectrum, V ne ; and Woo, et al. [1976b], have shown that if the 
three-dimensional wave number spectrum 4> of the electron density fluctua- 
tions has a power-law dependence with spectral index p ($> ne “ K~^ where K 
is the wave number), the frequency spectrum of the phase-difference 
scintillations and the one-dimensional frequency spectrum V ne (measured 
by in situ spacecraft) are both power-law with spectral indices p-1 and 
p-2 respectively. Specifically V Re is related to by 
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4.39 x 10 
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(18) 


where K ■ 2v/\ is the wave number, A is the wavelength, R is the closest solar 
distance of the line-of-sight path in units of R , v is the solar wind 
velocity transverse to the line-of-sight path at® R, f is the fluctuation 
frequency, c is the speed of light, F is the Gamma function, and all units 
are SI [Woo and Armstrong, 1979]. 
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FIGURE 11 

Frequency spectrum, W$d, for phase difference fluctuations shown in Figure 
10. The dashed line represents the power law f 8 “^’ 6 . The vertical bar 
defines the 90% confidence interval. The number of points used in the 
running point average m ■ 15 (after Woo et al. 1976b). 
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FIGURE 12 

Typical power spectra, W$<j, of the Viking phase-difference scintillations 
observed over the entire solar elongation (E) range. The spectra for 
R 5 30.5 RQare of open- loop data while the rest are of closed-loop data 
(after Woo and Armstrong, 1979). 
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This relation, (18), baaed on the Rytov geometric-optica aolution, 
appiiea for fluctuation frequencies leaa than the Freanel frequency, which 
la the order of 1 Hz for weak acintillationa when R > 20 R,. It allows the 
phaae spectra of Figure 12 to bo related to the spectra of electron density 
fluctuations. These are shown in Figure 13 for a large number of solar 
impact parameters. Note that the exponent of the fluctuation frequency 
exponent is now p-2 as compared to p-1 for the phaae spectra of Figure 12. 
Hence the best fit to the spectra of Figure 12 is 

W..(f) - 5.6 f" 2,65 R“ 2,45 2 (19) 

$d s Hz 

where R is in solar radii (20 < R < 100) and f is in Hz. 

3.2 Intensity Scintillations 

Interplanetary scintillations (IPS) of natural sources beyond our 
solar system have been studied extensively from networks cf observing sites 
on Earth's surface [Coles 1978, Coles and Harmon 1978]. The degree of 
scintillation is characterized by m, the Scintillation index, defined as 
the rms fluctuation of Intensity about the mean relative to the mean. 

Figure 14^^hows how it varies with the strength of turbulence and wave- 
length, because turbulence Increases as p decreases. Thus, m Increases 
with turbulence if the turbulence is weak, reaching a maximum < 1 at 

some . Increasing turbulence beyond this point produces a sharp 
decrease in scintillation, which is due to the angular width of the source. 
The more nearly the source approximates a point source, the cleser mm A ft 
approaches unity [Coles, Rlckett and Rutusey, 1974]. Figure 14b illustrates 
the behavior of the scintillation m f qjS. a point source (Hellos spacecraft) 
where it is seen that m becomes assymjNlitic to unity for R < 10 
[Woo, 1979]. 

For spacecraft transmission, the behavior or extraterrestrial scintilla- 
tion described in Report 263-4, and by Crane [1977], who has pointed out that 
the Interpretation of the spectral data in terms of signal level probability 
distributions has proved difficult. 

3.3 Spectral Broadening 

Physically, spectral broadening arises because of Doppler shifting 
of the rf signal as it is scattered by moving irregularities. Spectral 
broadening observations are similar to intensity scintillations In that 
they respond to^^nd furnish intc matlon on, only small-scale fluctuations. 
With spacecraft tRgnals, the closest distance to the Sun that can be use- 
fully probed by both spectral broadening and intensity scintillations is 
set by signal-to-roise ratio considerations and these are mere constraining 
for intensity scintillations [Woo and Armstrong, 1979]. In the case of 
Viking, spectral broadening was measured as close as 1.8 R # at 2.3 GHz and 
vt 1.4 R # at 8.4 GHz. Examples of spectral broadening for one-way traversal 


•i 



34 



FREQUENCY (H«) 


FIGURE 13 

Typical one-dimensional electron density spectra v ne deduced from Viking 
phase-difference scintillation spectra. Data marked A are in situ measure- 
ments from Pioneer 6 by Intriligator and Wolfe [1970], the curve marked B 
is from in situ measurements of 0G0 5 by Neugebauer [1975]. (After Woo and 
Armstrong, 1979). 




73.8 MHZ 


P (A. U. ) = SINE (ELONGATION) 


1977-1978 

HELIOS 


FIGURE 14 

Behavior of the intensity scintillation index, m, for natural (extended) 
and spacecraft (point) radio sources. (a) The variation of scintillation 
index with distance from the sun for natural sources at three radio fre- 
quencies. (After Coles et al. 1974). (b) Intensity scintillation measure 

ments of 2.3 GHz signal from Helios spacecraft received at Goldstone, 
California plotted as a function of nearest approach of the propagation 
path to the Sun (after Woo, 1979). 







at 2.3 GHz are shown in Figure 15. The spectral broadening bandwidth B 
is defined as in Goldstein [ 1 969 J : 



P(f) df 



P(f) df 


( 20 ) 


where P(f) is the spectrogram of the received signal. The radial variation 
of B with heliocentric distance at 2.3 GHz is given by Woo [ 1 978] . Spectral 
broadening has a serious effect on performance on paths with small impact 
parameters. 

3.4 Angular Broadening 

Decrease in signal level at very close approach to the Sun was observed 
both for Pioneer 6 [Goldstein, 1 969 J and for Helios [Woo, 1978]. As was 
discussed in Section 2.1, the collision frequency in the plasma is still 
too low to account for this decrease through absorption, so it is attributed 
to the angular broadening of the signal from the spacecraft due to turbulence 
to the extent that it is larger than the beamwidth of the receiving antenna 
[Woo, 1978]. In the case examined, it amounted to 4.5 dB and was clearly 
a significant consideration. 

3.5 Coherence Bandwidth 


There appears to be few If any experimental result" on phase and 
amplitude coherence of two monochromatic signals at different radio 
frequencies, but Woo [1975] has calculated the coherence functions for 
the X- and S-band cases. In these cases, coherence is maximum for signal 
fluctuation components of around 1 Hz. 


3.6 Frequency Dependence 


The frequency dependence for the homogeneous plasma parameters at 
microwave frequencies is determined by the simple asymptotic expression 
(f >> f N , f >•> f H ): 




( 21 ) 


These values are summarized in Table III in the main text, and, with the 
exception of dispersion (f”3), all vary with radio (carrier) frequency as 
f“2. 

The effects of Inhomogeneities in the plasma present a more compli- 
cated picture as their frequency dependence is a function of the turbulence 
parameter, p. These dependencies for the rms behavior of the various 
parameters are summarized in Table V. While the values listed are basically 
theoretical, the values of p are derived from the measurements in Figure 13 
where the exponent shown for v (f) is p-2. The value p = 11/3 is that 
correspond ing to the Kolmogorov spectrum. A striking feature of Table V is 
that only for Intensity scintillation does the expone... decrease in going 
from weak to strong scintillation. 
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TABLE V 


Frequency dependence of plasma turbulence effects in terms of the 
exponent y (effect ^fy where f is the carrier frequency) 





Exponent y 


Effect 

Parameter 

Relation 

Large Elongation 
(R * 20 Rg,) 
p - 11/3 

Small Elongation 
(R < 20 R e ) 
p ■ 3 

Intensity ... 

scintillation^ 

m 

. £_+. 1 
4 

-1.42 

-1.2 

Phase 

scintillation'' 

% 

-1.0 

-1.0 

-1.0 

Spectral 

broadening 

B 

2 

2 " P 

-1.2 

-2.0 

Angular 
broaden ing 

e 

P 

2 - p 

-2.2 

-3.0 


^Woo 1978, Woo et al., 1976a. 

(2) Woo 1978. 

^^Woo and Armstrong 1978, Woo et al., 1976a. 
(4) Woo et al., 1977. 
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SECTION III 


MODELS OF EXTRATERRESTRIAL PLASMAS 


1. INTRODUCTION 

It is clear from the discussion in Section II that three characteristics 
of extraterrestrial plasmas that affect microwave signals are electron 
density N, the one-dimensional spectrum of rms electron density fluctuations 
v ne’ an< * velocit y the pl asma motion (ordered and thermal). The 
effects of the solar corona and solar wind tend to be common to all space 
missions and these two topics are therefore treated in some detail in 
Section III. 2 and in Table VI. Space does not permit similar detail for each 
planet, so the principal features are tabulated in Table VII and references 
to some of the extensive literature on the subject are provided. 

It is not efficient to completely separate the discussion of propagation 
effects and of the propagation medium. Therefore, one of the three important 
characteristics of the propagation medium mentioned above has already been 
covered in considerable detail. This is the one-dimensional spectrum of 
rms electron density fluctuations (v ) that is treated in Section Ii«3. 
Hence, the part on inhomogeneities in this Section is very brief. 

2. THE SOLAR CORONA AND SOLAR WIND 

The characteristic of the Sun and terminology associated with this 
discussion of the Sun are summarized in Table VI. To the extent possible, 
the numerical values have been made consistent with those given by 
Allen [1976]. 


2.1 Plasma Density 


Berman [ 1 979] has recently reviewed some 16 determinations by 11 authors 
of the solar wind (or extended corona) term of the Allen-Baumbach relation 
[Allen, 1976] for the electron density out to 1 AU. The model now 
recommended by Berman is: 


N (R) 


2.21 x 10 


14 


1.55 x 10 

„2 . 3 


12 


electrons m 


-3 


( 22 ) 


In this model, the inner corona is modelled by the first term, which 
gives a value of N = 1.25 x 10 l4 m -3 at R = 1.1 R 0 , which agrees with 
Allen's [1976] model and becomes equal to the second term at 3.8 R @ . (For 
electron densities at lower solar altitudes, tabulated values should be 
used.) It yields a value at 1 AU (214.94 R 0 ) of N = 6.70 x 10^ electrons 
m -3 which is in good general agreement with observations, Allen [1976] 
points out that the electron density in the solar wind at 1 AU varies 
inversely with velocity and reaches a maximum of about 8 x 10^ m“ 3 . 
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TABLE VI 

Characteristics and terminology of the sun 


Heliographic equator: 

Rad ius , R s : 

Rotation period: 
(sidereal) 

Photosphere: 

Chromosphere: 


Corona: 


Solar wind 


The plane perpendicular to the Sun’s rotation axis 
that divides its mass in two; the inclination of the 
solar equator to the ecliptic is 7°15’ 

696,997 km 

25 days at equator 

28 days at 40° heliographic latitude 
>28 days nearer the poles 

Visible disk of the Sun; it has a physical temperature 
of 6000 K at its base, and a minimum of 4150 K at 
500 km above this base; the disk of the Sun subtends 
an angle at the Earth, which varies between 31 ’31" of 
arc and 32 ’35" of arc 

Spherical shell of 2000 km thickness actually made 
up of about half a million small spikes called 
spicules with lifetimes of 5 to 15 minutes; the 
chromosphere appears to have a temperature of about 
8000 K; all measures of solar radius, height in the 
corona, etc. , are made from the base of the 
chromosphere [Allen, 1976] 

Faint white halo around the Sun, which becomes 
visible during an eclipse; it is at an effective 
temperature of 2 to 6 million K; the term also refers 
to the actual plasma of electrons and ions in the 
upper solar atmosphere; the base of the corona may 
be taken as the transition region at 1.003 R@; the 
outer boundary is not clearly defined but is commonly 
taken at 6 R e [Yakovlev, 1974] 

The expansion of the corona into space; the corona 
contains areas that are particularly cool and quiet, 
known as coronal holes: many astronomers feel that 
the solar wind emanates from these holes; the solar 
wind originates inside the corona and is terminated 
at a few tens of AU by an interaction with the inter- 
stellar magnetic field [Brandt, 1970]; the solar wind 
velocity near the Earth is ^450 km/s and the average 
density at 1 AU is 5 to 10 x 10 6 m - -*; the magnetic 
field associated with the solar wind is initially 
radial but may be directed either inwards or outwards; 
a region of like-directed field is known as a sector 


TABLE VI (Continued) 


Coronal holes: Regions of diverging magnetic field of the corona, 

which appear dark when viewed in soft X-ray images, 
in white light comographs, and in eclipse photo- 
graphs of the solar corona; coronal holes are 
alwnys sources of high-speed solar wind streams; 
because coronal holes- tend to last longer than one 
solar rotation, there is a 27-day recurrence period 
for high-speed streams at the Earth, producing a 
natural explanation for the well-known recurrence 
period in geomagnetic activity [Bartels, 1934], and 
the associated source at the Sun, which has up until 
this discovery been known as an "M-Region" 

[Suess, 1979]. 
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The astronomical unit (the mean distance from the Earth to Sun) is 1.496 


Expression (22) is appropriate to low heliographic latitudes (±15°). 
Allen [1976] and Saito, et al., [1977] both find a decrease in electron 
density by a factor of 4 to 5 between the equator and pole for the inner 
corona. Tyler, et al., [1977] propose a model including a simple expression 
to account for variation with heliographs latitude (6), which has been 
modified slightly here to fit the above observations: 



1 2 \ 1/2 
6 + 25 sin 0 / 


(23) 


This expression may be used in connection with (22) if desired. As N(0) ■ 1 
no adjustment of (22) is necessary. Its appropriateness is questionable 
for R > 5 R e . 

It is now believed that the solar wind and, hence most interplanetary 
plasma, originates in coronal holes. A survey of models of coronal hole 
flows has recently been published by Suess [1979], while Holzer [1978] has 
written a review of the solar wind and related phenomena, which has been 
cast in the context of coronal holes. A general review of coronal holes 
and high-speed wind streams has been published by Zirker [l 977] . 

While substantial variation of plasma density with solar cycle is known 
to exist, these variations have yet to be sufficiently understood to be 
incorporated into models. For the time being, it is suggested that 
expression (22) be used for electron densities near the plane of the 
ecliptic, supplemented by expression (23) for higher heliographic latitudes 
for R < 5 R e , and that solar cycle variations be ignored. 


2.2 Inhomogeneities 


Fluctuations in radial electron density have been characterized in 
terms of the one-dimensional frequency spectrum V ne (measured in situ by 
spacecraft) by Woo and Armstrong [1979]. Their determination is: 

V (£) - 5.2 x I0 10 If 3 - 45 f-'- 65 V-' («l«tr.n./c.3)i 
ne s Hz 


and, if the solar wind speed v is assumed to be 3 x 10^ m/s and constant: 


V 

ne 


(f) 


1 . 5 x 


10 ' 


r-3.45 


3 

,-1.65 electrons/cm 
s Hz 


(25) 


f is the fluctuation frequency in 
s 


where R is in solar radii (20 R 215), 

Hz, and v is the solar wind speed in m/s. 

For values of R < 20 R , one can interpolate values of the V 
Figure 13. * 


ne 


from 


2.3 Velocity and Temperature 

The velocity vector follows the spiral defined by the rotation of the 
Sun and the frozen-in magnetic field [Brandt, 1970]. A recent scintillation 


analysis indicates that the average radial component of the solar wind 
velocity for distances less than 1 AU from the sun is given by 


V (r) - f" — Zj — km/second (26) 

r [r + 3.35 x 10 7 r * ,0 J 

where r is the radial distance in AU. 

There is a direct relationship between velocity and geomagnetic activity 
such that a velocity of 400 km/s is associated with very quiet and v - 700 
km/s with very disturbed geomagnetic conditions [Allen, 1976]. 

3. PLANETARY IONOSPHERES AND MAGNETOSPHERES 

This section discusses the four inner planets (Mercury, Venus, Earth, 
and Mars) and the nearest (and largest) of the outer planets, Jupiter. At 
present these are the planets for which data from spacecraft visitations 
are available. Summary information is found in Table VII. The atmospheric 
information is drawn from Sagan [1975], and Pasachoff [1977]. Information 
on planetary ionospheres is taken from the sources referenced in the texts 
on the various planets in this section and its associated figures. 

3.1 Mercury 

Mercury, the innermost planet, has one of the least hospitable surface 
environments, due to its close proximity to the* Sun. The maximum angle from 
the Sun at which Mercury can be viewed from the Earth is only 28° so that it 
is never seen optically against a really dark sky. Until recently astronomers 
believed the rotation period of Mercury to be the same as its period of 
revolution about the Sun, namely 88 days. This would have meant that one 
half of the planet was always in sunlight, the other half always in dark. 
However, radar measurements from the Arecibo Gbservatory indicate the 
rotation period to be 59 days. Thus the spin-orbit coupling, instead of 
exhibiting a ratio of 1:1 has a ratio of 2:3. Further information may be 
found, for example, in Murray [1975] and Pasachoff [1977]. 

The existence of a magnetic field on Mercury has been established and 
electrons and traces of atmospheric gases observed. The planet exhibits a 
bow shock and magnetopause similar to that of the Earth, but both are 
severely compressed in the direction of the Sun. 

? , 2 Venus 


Venus is covered by clouds from about 30 km above the surface to an 
altitude of 80 km. The visible surface of these clouds at around 60 km 
exhibit a temperature of around 240 K at a pressure of around 200 mb. From 
this level, the temperature increases monotonically to a value of about 750 K 
on the solid surface, at which point the pressure has risen to 90 bars. The 
clouds of Venus give the planet an apparent rotation period of about four 
days [Young and Young 1975], while radar measurements of the solid surface 
put the rotation period at 243 days. The difference can be explained by 
postulating wind velocities in the cloud layer of up to 100 meters per second. 
Good discussions are found in Young and Young [1975J, Pasachoff [1977], and 
Colin, et al., [1 979]. 
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1 1 -3 

Electron densities of 5 x 10 el m have been measured at altitudes 
of 150-180 km above the solid surface [Knudsen, et al., 1979]. The ioniza- 
tion is found to stay about 10 ^ e i m - ^ for ionopause heights of 1500 km, 
and then to drop sharply above the ior.opause. The principal positive ion 
above 180 km is 0 + and below 180 km it is 02 + [Taylor, et al., 1979]. Some 
of the measured data on Venus is given in Figure 16. The maximum electron 
density shown corresponds to a plasma frequency of 6.3 MHz [Kliore, et al. 
1979' 

1979 


3 : 


The magnetic field has been estimated to be very weak [Russell, 


3.3 The Earth 


The ionosphere of the Earth up to an altitude of 2000 km is described 
in detail in Report 725, while the effects of Earth-space propagation out 
to the geostationary orbit at 6.6 Earth radii measured from the center of 
the Earth (6.6 IL_) are given in Report 263-4. In this section, the termi- 
nology necessaryto distinguish ionization connected to the Earth from that 


of the solar wind is given. 

The Earth's plasma environment Is shown in Figure 17 [Van Allen, 1975]. 
This environment is strongly influenced by the Earth's magnetic field repre- 
sented by a dipole moment of 8 x 10“ gauss cm“^ inclined 11“ t 
of rotation. 


to the axis 


3.3.1 Plasmasphere 

O i n 

This region of relativity high electron densities (10’ to 10 electrons 
m”^) at low intermediate geomagnetic latitudes is contained along closed 
lines of geomagnetic force. In the equatorial plane, the plasma pause occurs 
between 5 and 6 Rg* and defines an asymmetrical torus outside of which the 
electron density is 10 to 100 times less than Inside. The location of the 
plasmapause is strongly influenced by magnetic »ctivity. 

3.3.2 Magnetosphere 

In this region, the Earth's plasma is controlled by the geomagnetic 
field. The boundary between the ionosphere and magnetosphere is taken by 
some authors to lie at altitudes of 500 to 2000 km above the surface, while 
other authors choose not to distinguish the regions. The magnetopause 
is the outer boundary of the magnetosphere and is seen to be compressed on 
the day side of the Earth and extended on the night side. 

3.3.3 Bow shock 

The bow shock separates the undisturbed supersonic solar wind from the 
bow-shaped shock wave that lies between it and the magnetosphere. The 
electron density in the solar wind impinging on the Earth's bow shock ir 
approximately 5 to 10 x 10^ el m - ^. As in all plasmas, electric neutrality 
is maintained so there are an equal number of positive ions. 
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FIGURE 16 

The ionosphere of Venus from Pioneer Venus measurements. (a) Electron den- 
sity measurements acquired by the orbiter electron density probe on 1978 
day 356, showing the bow shock, A, A', and the ionosphere, B. (After Brace 
et al. 1979); (b) Ion composition during entry near dawn local time 
(Taylor et al. 1979); (c) A series of electron density measurements of 

orbits 1, 4, 9, 13, 17 and 21. The profiles are staggered by one decade, 
and the 1(P em - ^ level is indicated for each profile. The solar zenith 
angle (x) for each orbit is indicated on the upper horizontal axis. 

(Kliore et al. 1979). 
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3 . 4 The Moon 




The Moon is the most massive planetary satellite in the solar system 
in relation to its planetary partner. It is without atmosphere or ionosphere 
The Moon's mean distance from the Earth, about 380,000 km (or about 60 R^) , 
and its orbital inclination are such that the Moon moves through the tail 
of the Earth's magnetosphere once a month at full Moon. 

3.5 Mars 

Mcrs is a small planet only about half the diameter of Earth. The 
Martian c.'-mosphere is about one percent that of the Earth (around 6-mb 
pressure at the surface), and is 90-percent carbon dioxide. The surface 
temperature is below freezing (-23°C) and there is, in fact, little water 
vapor in the atmosphere. The Martian magnetic field is very weak (see 
Table VII). 

Shown in Figure 18 are two sets of measurements (Vikings 1 and 2) of 
the ionosphere of Mars [Hanson, et al., 1977]. The peak electron density 
(5 x 10^ m“3) corresponds to a plasma frequency at 2.85 MHz and is 
significantly less than that of the Earth's F-region. 

3.6 Jupiter 

Jupiter, the largest of the planets, dominates the Sun's planetary 
system at a distance of 5 AU. It revolves around the Sun once every 
11.9 years. Its volume is 1400 times that of Earth, but its mass is only 
318 times as great (its density is only 1.3 g/cm^) . It rotates very 
rapidly (once every 10 hours). Jupiter's atmosphere is 85-percent hydrogen 
molecules and 15-percent helium, plus traces of methane and ammonia. There 
are at least 14 moons, four of them (discovered by Galileo) of a size 
comparable to the planet Mercury: Io, Europa, Ganymede, and Callisto 

[e.g. , Wolfe, 1 975 J . 

That Jupiter has a strong magnetic field and radiation belts similar 
to the Earth was deduced from noise bursts emanating from Jupiter and 
measured on Earth by Burke and Franklin [ 1 955] at 22 MHz. Since then, 

Jupiter has receivtd a lot of attention from radio astronomers primarily 
at decametric wavelengths [see Kraus, 1966], Recent measurements of the 
Jovian ionosphere are shown in Figure 19. 

3.7 Other 

The planets Saturn, Uranus, Neptune, and Pluto are not covered in this 
report, other than the sketchy statistics included in Table VII. The Pioneer 
Venus spacecraft encountered Saturn in September of 1979, and new information 
on that planet is expected shortly. Scientific American [ 1 975] offers 
general information on the outer planets. 

Comets offer an interesting future area for spaceflight study. The 
tail of a comet first becomes visible about 1.7 AU from tr.e Sun [Allen, 1976] 
and frequently separates into two distinct components: an ion tail that is 
always pointed away from the sun, and a dust tail that curves along the 
comet's parabolic trajectory around the Sun [Pasachoff, 1977]. 
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APPENDIX 


AUTHORIZING TEXTS OF CCIR STUDY GROUP 6 


At the CCIR XIV Plenary Assembly, Kyoto, 1978, the terms of reference for 
Study Group 6 (Propagation in Ionized media) were expanded, a Question was 
broadened, and a Study Programme was approved to Include propagation through 
extraterrestrial ionized media. In CCIR usage the Question and Study Programme 
are defined as follows: 

Question : A Statement of a technical or operational problem, to which an answer 

is required. 

Study Programme : Text describing the work to be carried out on a technical or 

operational problem constituting the subject of a Question. 

The specific wording approved at the 1978 Plenary Assembly is: 

Terms of reference (for Study Group 6, Propagation in ionized media) to study: 

1. the propagation of radio waves through the ionosphere, and through ionized 
regions beyond the ionosphere; 

2. the characteristics of radio noise; 

with the object of improving telecommunications. 


QUESTION 25/6 
IONOSPHERIC PROPERTIES 

The CCIR, (1978) 

CONSIDERING 

that ionized media affect the propagation of radio waves, 

UNANIMOUSLY DECIDES that the following question should be studied: 

what are the properties of the ionosphere and other ionized media which 
affect propagation in ways that are important to radio systems? 

Note. - See Reports 250-4, 262-4, 340-3, 725, 726 and Recommendations 373-4, 
434-3. 
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STUDY PROGRAMME 25E/6 


RADIO PROPAGATION THROUGH THE SOLAR 
AND OTHER EXTRA-TERRESTRIAL IONIZED MEDIA 


The CCIR, (1978) 

CONSIDERING 

(a) that a simple model of the ionized solar atmosphere and the ionospheres of 
the planets would provide a basis for propagation studies needed to support the 
planning of space flights near the sun and the planets; 

(b) that simple models of the extra-terrestrial plasma would provide a similar 
basis for deep-space flights, 

UNANIMOUSLY DECIDES that the following studies should be carried out: 

1. determination of the characteristics of the ionized solar atmosphere, the 
ionospheres of the planets and their associated effects on radio propagation; 

2. determination of the characteristics and effects of the extra-terrestrial 
plasma pertinent to radio propagation over interplanetary to intergalactic distances. 
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